from scipy.interpolate import interp1d

# Define data points
x = [1, 2, 3, 4, 5]
y = [2, 3, 5, 8, 12]

# Create an interpolation function
f = interp1d(x, y, kind='cubic')

# Evaluate at new points
new_x = [1.5, 2.5, 3.5]
print("Interpolated values:", f(new_x))